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Abstract 

We point out that the recent Belle measurements of the B d j s — > J/ijnf'' decays imply large 
pseudoscalar glueball contents in the rf 1 ' meson. These decays are studied in the perturbative 
QCD (PQCD) approach, considering the n-n'-G mixing, where G represents the pseudoscalar 
glueball. It is shown that the PQCD predictions for the B d / S — > J/iprj^ branching ratios agree 
well with the data for the mixing angle <pc ~ 30° between the flavor-singlet state and the pure 
pseudoscalar glueball. Extending the formalism to the n-n'-G-ric tetramixing, the abnormally 
large observed B d — > Kn' branching ratios are also explained. The proposed mixing formalism 
is applicable to other heavy meson decays into rfi^ mesons, and could be tested by future LHCb 
and Super- B factory data. 
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Recently, the Belle Collaboration reported their new measurements of the B d — >■ J/ Wr/ r> 
decays and the first observations of the B s — > J/ip-qW decays j2[ with the branching 
ratios 

BR(B d J/^) Exp = 12.3: 
BR(B d -> J/HWv < 7.4 
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BR(5 S -> J/ipvhxp = 5.10 ± 0.50(stai) ± 0.25(syst. 
BR(5 S -+ J/iprfhxp = 3 - 71 ± 0.61(stat) ± 0.18(syst. 
and the relation of the branching ratios |2[ 
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The updated result in Eq. (JTj), as stated in [l|, is consistent with and supersedes the 
previous one, BR(B d — >■ J/iprj) = 9.5^\^ x IO -6 , in 2007 [3]. The accuracy of the above 
data is expected to be improved rapidly along with operation of LHCb experiments. 

The B d /s —> J/iprj^ decays were first investigated by Deandrea et al. j3] under the naive 
factorization assumption without considering the 77-77' mixing. In Ref . [5| , Skands analyzed 
the B d / S — > J/tprj decays via SU(3) relations to the "golden channel" B d — > J/ipK° 
including the 77-77' mixing. The predictions BR(B d — > J/iprj) = 11.0 x 10~ 6 and BR(5 S — > 
J/ijyq) = 5.0 x 10" 4 {BR(B d J/ijyq) = 15.0 x 10~ 6 and BR(B S -»■ J/^jr]) = 3.3 x 10~ 4 ), 
corresponding to 6 = —10° (6 = —20°), were obtained [5|, where 9 is the angle between 
the flavor- singlet state 771 and the flavor-octet state r)&. Note that the B d / S — >■ J/iprj' 
branching ratios were not calculated in If they were, the outcome of BR(B S — > J/ijnf) 
would be much larger than in Eq. (j3J). This must be the case in the conventional r\-v( 
mixing formalism as elaborated below. 
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FIG. 1: (Color online) Leading quark-level diagrams for the B d u — > J/iprj^ decays. 



The 77W mesons are produced via their nonstrange components in the B d meson decays 
at leading level, while they are produced via their strange components in the B s meson 
decays, as easily seen in Fig. [TJ Under the conventional mixing, the physical states rj and 
7/ are related to the quark-flavor states rj q = (uu + dd) / \/2 and r] s = ss through a mixing 
matrix with an angle (f) j6[, 
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Various experimental and theoretical constraints have yielded <fi < 45° (for a recent de- 
tailed discussion, see Ref. 0), for which one observes (i[ 

s w < 1 , (7) 

R s = ^7— — r- « cot 2 > 1 . 8 

BR(B S ->■ J/^77) ^ v ; 

It is obvious that only the first relation meets the data, and the second one does not as 
confronted with Eq. (jSJ). As speculated in Q, the small measured ratio Rf xp indicates 
additional flavor-singlet components in the r\ meson other than the uu, dd and ss pairs, 
or violation of the 77-77' mixing scheme. 

Inspired by the Belle measurements, we analyze the B^/s J/ipr]^ decays in the 
more complete ryn'-G mixing formalism, G denoting a pseudoscalar glueball. Adopt the 
construction in [9(, 

'\v)\ ( \m) \ 

\rf) =U^e)U 1 {4> G ) It?!) , (9) 
. \°) J \ \9) ) 

where the matrices 

10 \ ( cos9 -sinO 

Ui(<p G ) = I cos <pQ sin0 G , U 3 (9) = sinfl cosfl | . (10) 
— sin (pa cos 4>g ) \ 1 

represent rotations around the axis along the 7/g meson and the unmixed glueball g, 
respectively. Equation (Q is based on the assumption that r/ 8 does not mix with the 
glueball, under which two mixing angles 9 and <fi G are sufficient. 

The octet and singlet states are related to the flavor states through 



fill 



where 9{ is the ideal mixing angle with cos 6^ = yT/3, i.e., 9\ = 54.7°. The flavor states 
are then transformed into the physical states via the matrix 

[/(</>, <j> G ) = U^U^gPM , 

cos cf) + sin 9 sin 9{ A G — sin + sin 9 cos 9i A G — sin # sin 0c* \ 
sin — cos 6* sin #j Ac cos — cos 9 cos 9i A G cos 6* sin cj) G , (12) 
— sin 0i sin (j) G — cos 9{ sin cj) G cos 4> G J 

with the angle <fi = 9 + $i and the abbreviation A^ = 1 — cos<f) G . U has been written in 
the form, that approaches the mixing matrix in Eq. (J6]) in the <p G — > limit. We mention 
that a different mixing matrix has been constructed by diagonalizing a mass matrix for 
pseudoscalar bound states, which was derived from an effective QCD Lagrangian (Io| . 

It has been verified that the contribution from the gluonic distribution amplitudes in 
the 77W meson is negligible for B meson transition form factors [11] . Hence, the 77 and 77' 
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mesons are still produced via the nonstrange (strange) component in the B d (B s ) meson 
decays under the rj-rj'-G mixing. Equation f[T2"j) then simply modifies Eqs. (j7|) and OH]) into 
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Th 



sin — cos sin 0j Aq 
cos + sin 9 sin 6*j A G 
cos — cos # cos 6i Aq 
— sin + sin 6 1 cos 9i Aq 
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respectively. Since the angle 9 is negative, both the numerator and the denominator 
in Eq. (1131) decrease, and the ratio Rj h could remain smaller than unity. On the other 
hand, the numerator in Eq. ( |14p decreases, while the magnitude of denominator increases. 
Therefore, the ratio Rj h might drop from above unity to below unity for a sufficiently 
large angle 4>g- 

We compute the B d / S J/ipr/^ decays explicitly in the rj-r/'-G mixing formalism, 
employing the perturbative QCD (PQCD) approach [12|] at next-to-leading order (NLO) 
of the strong coupling constant. Replacing the kinematic variables and distribution am- 
plitudes of the K° meson by those of the ru s ) meson, we derive the Bd — > J/ipr/^ decay 
amplitudes from the B d — > J/tpK ones [l3j j. Further substituting B s for B d , we arrive 
at the B s — > J/ipr)('> decay amplitudes. Throughout this work it is assumed that there 
are no final state interactions and no isospin violation. The input parameters such as 
the QCD scale (in units of GeV), masses (GeV), decay constants (GeV), and B meson 



lifetime (ps) are set to [12|, [l 
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(15) 



Following the observation in [15(, we have chosen a value for m slightly higher than the 

1.3 ~ 1.4 GeV. The best fit on the mixing angle 9 from 



chiral scale for the pion, 
experiment data converged in the range of [—11°, —17°] 16]. Our choice = 43.7° in 
Eq. (j!5p is equivalent to 9 = — 9i = —11°, close to 9 = —10° in [5], which corresponds 
to a sizable gluonic admixture. It is also consistent with = (44 17) determined from 
the B d — > J/^prj/n data [l7|. The angle 0g, as summarized in [9j, varies in a wide 
range 10° < 6n < 30°. We have taken the central value of 0g = 33° ± 13° recently 
extracted in |l7J, |l8| . For the Cabibbo-Kobayashi-Maskawa matrix elements, we employ 
the Wolfenstein parametrization with the updated parameters A = 0.832, A = 0.2246, 
p = 0.130, and fj = 0.350 Ei . 



The CP-averaged branching ratios for the B d / S J/ipr/^ decays in the standard model 
then read as, 



BR(B d ->■ J/ijf]) 
BR(B d -> J/H) 
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with the ratios Rj h ~ 0.58 and Rj h ~ 0.72. We have kept the 77' meson mass in the 
phase-space factor for the B d / S — >■ J/iprj' decay rates, and neglected masses of other light 
mesons. The theoretical uncertainties arise from the variation of the shape parameter 



tu B = 0.40 ± 0.04 (tu Bs = 0.50 ± 0.05) GeV for the B d (B s ) meson wave function |19|, |20 
of the J/ij) meson decay constant fj/^ = 0.405 ± 0.014 GeV [2lf] , and of the Gegenbauer 
coefficient a<i = 0.44 ± 0.22 for the leading-twist rj q ^ distribution amplitude. Obviously, 
the theoretical branching ratios are in good agreement with the existing data and upper 
bound after considering the rj-rf-G mixing. 

A remark is in order. The B s — > J/iprj^ decays have been proposed [22} to explore the 
i]-i]'-G mixing with KLOE's parametrization for the mixing matrix (23[, namely, assuming 
the absence of the glueball content in the rj meson. It was observed that the angles ~ 40° 
and <f>o re 20° determined by KLOE lead to R s ~ 1, consistent with the previous Belle 
data [24| . Hence, the updated data R s < 1 indeed imply larger glueball components in the 
rfi'> meson as explained before. In principle, R s < 1 can be achieved in both KLOE and 



our parametrizations by tuning the angles. Our parametrization is close to that in [17J, [18 



with a nonvanishing glueball component in the 77 meson, so we chose the larger <j)Q = 33° 



extracted in [17j, [18(, and obtained R s < 1 naturally. Taking <p G = 22° [23J, the central 



values of the branching ratios become 

BR(B d J/iprj) = 11.7 x 10~ 6 , BR(B d J/ipr]') = 8.2 x 10" 6 , (20) 
BR(B S J/tjyq) = 5.00 x 10" 4 , BR(B S ->• J/ifrrf) = 4.28 x 10" 4 , (21) 

and the consistency with the data deteriorates. Compared to Eqs. f[T7[) and f[T§l) . we 
find that BR(B d / s — >• J/iprj) are less sensitive to (pc than BR(B d / s — > J/iprj'), a result 
attributed to the smaller glueball component in the 77 meson. Another remark is that as 
far as the B d — > J/ipr/^ decays are concerned, their branching ratios can be accommodated 
in the conventional 77-77' mixing by tuning the angle 0. Fitting Eq. ([7]) to our predictions 
in Eqs. and dUJ yields w 37.3°, close to = 39.3° ± 1.0° in (J. 

It is interesting to examine whether D, D s decays into rf'> mesons, such as D, D s —t 
i]^£ + i>, reveal the similar implication on the mixing mechanism 0,0, 26]. The ratios of 



the branching ratios of these semileptonic decays have been expressed as 

Rd = BR(D + -+ nl+v) ~ RDtaU ' (22) 
BR(D S -> n'l + v) 7 , . 

K ^ BR(Z)U VW <23> 

where the factors Rd ~ Rd collect the information on the D d / S — > i] q / s transition form 
factors and the corresponding phase space. The estimate of Rd depends on how to model 
the q 2 dependence of the form factor, q 2 being the lepton-pair invariant mass squared, and 
suffers theoretical uncertainty. Taking Rd ~ 0.28 j27|, it was found that the measured 



values Rf xp w 0.19 ± 0.05 < R D H and Rf xp = 0.37 ± 0.10 > R D Q exhibit a 
pattern in agreement with the conventional 77-77' mixing. However, we stress that the 
above observation is not in conflict with the f]-r]'-G mixing formalism advocated in this 
work, viewing the potential uncertainties in the estimate of Rd and in the assumption 
Rd ~ Rd- I n addition, it is not sure that the contributions from the D,D S transitions 
to pseudoscalar glueballs, which were referred as the weak annihilation process in |7J], are 
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negligible as in the Bd/ S meson decays. The inclusion of these channels would modify the 
analysis of the D, D s — > r]^t r v decays. 

At last, we have confirmed that the larger angle 0g = 33° generates the pseudoscalar 
glueball mass uig w 1.49 GeV, following the rj-rf-G mixing formalism in j^]. Namely, the 
postulation ^ that the the 77(1405) meson is a leading candidate for the pseudoscalar 
glueball is not altered. A larger <pG also makes an impact on the B — > Krf* branching 
ratios in the r)-r)'-G-r] c tetramixing formalism [15J| . Substituting (pa = 33° into the formula 
in Sec. Ill D of Q, BR{B° -»■ K°if) is enhanced from 50 x 10" 6 in NLO PQCD @ to 
(59. 7^4) x 10~ 6 . On the other hand, BR(_B — > K°rj) is insensitive to the variation of 
t>G, which remains as around 2 x 10~ 6 [15]. Therefore, the dramatically different data for 



can 



BR{B° -> K°rf) = (66.1 ± 3.1) x 10~ 6 and BR{B° -> K°rf) = (1.12+°;^) x 10~ 6 [30 
also be understood. 

In summary, most of studies in the literature on B^/s decays into r/^ mesons were 
performed under the conventional 77-77' mixing. The recent B^/s — * J/ipr/^ data provided 
a strong implication on the sizable pseudoscalar glueball contents in the ry>'> meson, which 
motivated our investigation in the 77-7/-G mixing formalism. We have verified this impli- 
cation by computing explicitly the B^/s J/ipV^ branching ratios in the NLO PQCD 
approach: the outcomes from a large angle 4>g ~ 30° were found to be well consistent with 
the current measurements and upper bounds. The abnormally large observed B — > Krj' 
branching ratios were also accommodated in the rj-rf-G-rjc tetramixing formalism with the 
same (j>a- Our work suggests that complete understanding of dynamics in 77 ^ -involved pro- 
cesses demands the rj-rj'-G mixing scheme. The resultant predictions for other B^/s — > rft' 
decays could be tested by future data of LHCb and/or Super- B factories. 
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